Ultrastrong coupling regime of cavity QED with phase biased flux qubits 
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We theoretically study a circuit QED architecture based on a superconducting flux qubit directly 
coupled to the center conductor of a coplanar waveguide transmission-line resonator. As already 
shown experimentally [Abdumalikov et al. Phys. Rev. B 78, 180502 (2008)], the strong coupling 
regime of cavity QED can readily be achieved by optimizing the local inductance of the resonator in 
the vicinity of the qubit. In addition to yielding stronger coupling with respect to other proposals 
for flux qubit based circuit QED, this approach leads to a qubit-resonator coupling strength g which 
does not scale as the area of the qubit but is proportional to the total inductance shared between 
the resonator and the qubit. Strong coupling can thus be attained while still minimizing sensitivity 
to flux noise. Finally, we show that by taking advantage of the the large kinetic inductance of a 
Josephson junction in the center conductor of the resonator can lead to coupling energies of several 
tens of percent of the resonator frequency, reaching the ultrastrong coupling regime of cavity QED 
where the rotating-wave approximation breaks down. This should allow an on-chip implementation 
of the E ® P Jahn- Teller model. 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Ta 



I. INTRODUCTION 

Combined with the large electric dipole moment of 
superconducting charge qubit, the large vacuum elec- 
tric field of microwave transmision-line resonators can 
be used to reach the strong coupling regime of cavity 
QED However, charge qubits suffer from charge fluc- 
tuations which leads to low coherence times. By working 
with a Cooper Pair Box qubit in a parameter regime 
where charge dispersion is small, the transmon qubit vM 
has led to significant improvement in coherence times [3| , 
in addition to larger qubit-field coupling strengths g. 
This is, however, done at the cost of lower anharmonicity, 
limiting gate speed. In Refs. it was suggested that 

an alternative approach to reaching the strong coupling 
regime with superconducting qubits is to inductively cou- 
ple flux qubits to the zero-point motion magnetic fleld of 
a transmission-line resonator. In this case, coupling in- 
creases with qubit loop area A, with A ~ 8 /im^ expected 
to be sufficient to reach coupling strengths of a few tens 
of MHz 0. While comfortably in the strong coupling 
limit, the predicted values are almost an order of mag- 
nitude lower than what can be obtained with transmon 
qubits Larger couplings can be obtained by increas- 
ing the qubit area, but only at the expense of increased 
sensitivity to flux noise. 

In this paper, we theoretically investigate an approach 
experimentally realized by Abdumalikov et al. [6[ where 
flux qubits are directly connected to the center conduc- 
tor of a coplanar waveguide transmission-line resonator. 
By changing the width of the center conductor to take 
advantage of the kinetic inductance, the phase bias of the 
qubit by the resonator is enhanced. We show how this 



approach leads to signiflcant qubit-resonator coupling, 
easily reaching the strong-coupling regime. Inserting a 
Josephson junction in the center conductor of the res- 
onator, much stronger couplings can be obtained, with 
g reaching several tens of percent of the resonator fre- 
quency. In this ultra-strong coupling regime, the ubiq- 
uitous rotating-wave approximation is expected to break 
down, leading to as of yet unexplored physics in cavity 
QED. In addition to the larger coupling, an advantage 
of this approach over that presented in Ref. [1, Q is that 
g does not scale with the qubit area. Moreover, with 
its multi-level structure, the flux qubit can be used in 
the A configuration Q opening the possibility to realize 
EIT in a cavity Q and a wealth of other quantum optics 
phenomena circuit QED. 

This paper is organized as follows. We start by finding 
the normal modes of an inhomogeneous transmission-line 
resonator. The case of a Josephson junction playing the 
role of the inhomogeneity is then discussed. Building on 
these results, we obtain the Hamiltonian for a flux qubit 
directly connected to the center conductor of the inho- 
mogeneous transmission-line and obtain expressions for 
the qubit-resonator coupling strength. Finally, numerical 
results for the coupling strength are presented. 



II. ENHANCED PHASE BIASING 

A schematic of the circuit we consider is shown in 
Fig. [T] A superconducting flux qubit is fabricated such 
that its loop is closed by the center conductor of a 
transmission-line resonator of length 21 (ranging from 
X — —i to Assuming a non-imiform resonator, the 



2 




Figure 1: (Color online) Schematics of a 3-junction flux qubit 
directly connected to the center conductor of an inhomoge- 
neous superconducting transmission-line resonator. 



Lagrangian density reads 

C°(x)^^(x,t) 
i-ti = 



1 fdip{x,t) 
2L"{x) \ dx 



(2.1) 



with %l){x,t) = j*_^dt'V{x,t'), C°{x) the position- 
dependent capacitance per unit length and L^(x) = 
-^geo(^) + -^kinetic (^) posltion-dcpendent inductance 
per unit length including both geometrical and kinetic 
contributions. 

The corresponding Euler-Lagrange equation of motion 



d 

dx 



1 d->p{x,t) 



L°{x) dx 



C"{x)4){x,t) (2.2) 



is solved by first decomposing ip{x, t) over (unitless) nor- 
mal modes u„{x)^ 



ip{x,t) = y^?/'„(t)M„(x). 



(2.3) 



Here ■(/)„ is the flux amplitude of eigcnmode n, of fre- 
quency uJn and eigenfunction Un{x) and is given by 

1 /■+^ 

V-nW = ^ C\x)^j{x,t)Un{x)dx, (2.4) 

where iV is a normalization constant. Assuming a large 
quality factor Q, the current at the two ends of the res- 
onator vanish; the eigenmodes must satisfy the boundary 
conditions d^u^^x = ±£) = 0. Spectral decomposition of 
the flux ip{x,t) in Eq. (|2.2p leads to a Sturm-Liouville 
differential equation of the form 



d_ 

dx 



1 dUnjx) 

L^{x) dx 



= -c.2c°(x)w„(x), 



(2.5) 



whose solutions {un{x),ujn} form an orthogonal basis. 
The eigenfunctions u„ (x) respect a weighted orthogonal- 
ity relation 

/ C°{x)Un{x)Ujn{x)dx = CrSnm, (2.6) 
J-t 

where the normalization constant is chosen to be the total 
capacitance of the transmission-line, Cr — J_f. C^(x)dx. 
By using the spectral decomposition (|2.3p in the La- 



- Si„= 5 m, ti„= 200 nm 
■ Si„= 50 nm, ti„= 200 nm 
• Si„= 50 nm, ti„= 50 nm 




Figure 2: (Color online) First mode (normalized to 1) for 
aluminum transmission-line resonators with the different con- 
striction parameters listed in Table H] As the central line is 
reduced in cross section, the slope of the flux field increases 
inside the constriction. The inset shows the geometry of the 
constriction in the center conductor of the resonator. 



grangian density (j2.ip and using the orthogonality re- 
lations (|2.6p along with the Sturm-Liouville differential 
equation (|2.5p . the total Lagrangian simplifies to a sum 
over eigenmodes: 



ECr 
~2 



I A -I 2/2 



(2.7) 



Defining the charge On = Cripn as the conjugate momen- 
tum to the flux ipn, the corresponding Hamiltonian is 



3_ 



2 ; 2 



(2., 



By quantifying and introducing the operators 



" V — 2 — 



(2.9) 



with [a„, aj„] = ^nm, we arrive at the standard form 



TYti = ^ tiwn (ajja„ + 1/2) 



(2.10) 



completing the mapping of the inhomogeneous resonator 
to a sum of harmonic oscillators. Unlike the homoge- 
neous case, the mode frequencies can be inharmonically 
distributed such that the equality u;„ = moo is not sat- 
isfied in general. Like most Sturm-Liouville problems, 
the eigenmodes Un{x) and eigenfrequencies w„ are found 
numerically by exact diagonalization Q. As discussed 
in Appendix [XI details of the transmission-line geometry 
are important in determining these quantities. 



3 



Table L Inductance ratio L^n/L^ut inside and outside the constriction, resonant frequency ijJi/2-k of the first mode, fiux gradient 
inside the constriction I^Vi {x) /9x|^^g and qubit-resonator flux coupling g^-^ at <^e.xt = $o/2 for different values of the resonator 
center conductor width Sin and thickness tin at the location of the constriction assuming a length of = 5 /^m of the shared part 
between the qubit and resonator. Qubit parameters are given in the text while resonator parameters are given in Appendix [XI 
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Material 


(nm) 


(nm) 


rO /rO 


(GHz) 


(10-6$o/Aini) 


(MHz) 


Al 


5000 


200 


1 


13.12 


12.97 


71.8 


Al 


50 


200 


3.4 


10.98 


35.36 


195.6 


Al 


50 


50 


4.1 


10.52 


40.22 


222.5 


Nb 


50 


50 


8.3 


8.62 


65.10 


360.2 



Figure [2] shows the first mode ui{x) for three different 
configurations of a constriction in the center conductor of 
the resonator, as detailed in Tabic [D As the constriction 
is made narrower, and thus the local inductance made 
larger, an abrupt change in ui{x) develops. A flux qubit 
connected on either side of the constriction, as illustrated 
in Fig. [31 will thus be strongly phase biased. As a result, 
an inhomogeneity in the resonator can increase the qubit- 
resonator coupling, making the strong coupling regime 
easier to reach. 

Qubit-resonator coupling can also benefit from the ki- 
netic inductance of superconducting materials with large 
London penetration depth Al, such as niobium where 
Al = 39 nm, when the dimensions of the cross-section of 
the central line reaches dimensions comparable with Al. 
Alternatively, a Josephson junction with large Joseph- 
son inductance Lj ^ ip^/Ej can replace the constriction 
shown in Fig. [T] to provide even stronger coupling. In 
Appendix [Bl wc show how in this case, the field ipi^i t) 
becomes discontinuous at the location of the junction and 
presents an important fiux difference across the junction. 
Wc also show how the Hamiltonian of the transmission- 
line plus junction can also be written in the standard 
form of Eq. (|2.10p [lo| . Consequences of this very large 
coupling will be discussed further below. 



III. QUBIT-RESONATOR HAMILTONIAN 

In this section, we obtain the qubit-resonator Hamil- 
tonian for the system of Fig. [1] focusing on the case of 
the inhomogeneous resonator. Derivation of the Hamil- 
tonian in the presence of a fourth junction instead of a 
constriction can be done following these lines and is dis- 
cussed in Appendix[Bl Figure [3| shows in more details the 
qubit connected to the center conductor of the resonator. 
Including the flux qubit, the Lagrangian reads 

^^^^^ + ^J^"^'- + ^°^['^'^-/^o]} , (3.1) 

where Cjk is the capacitance of junction k, Ej}. its 
Josephson energy and (pk the flux difference across it, 



lp(x,) 1|J(X2) 




Figure 3: (Color online) Closeup of the flux qubit fabricated 
at the location of the constriction. The qubit is attached at 
positions xi and X2 on the resonator. $cxt is an externally 
applied flux. A large Josephson junction inserted in the center 
conductor of the resonator between x\ and X2 can lead to 
stronger coupling. 



and LfQ = <I>o /2tt is the reduced flux quantum. Junctions 
1 and 3 are assumed to be equivalent, Cji = Cjs = Cj 
and Eji ~ Ejs = Ej, while junction 2 is such that 
Cj2 = aCj and Ej2 = aEj with a < 1 [H. The flux 
differences 4'i{2) depend explicitly on the resonators volt- 
age through "0(2^1(2))! with xi and X2 the positions of the 
resonator where the qubit loop is connected. Finally, the 
phase differences satisfy 

(^1 +(^2 -(^3 - V^tl = $oxt, (3.2) 

where "^ti = il'{x2) — ''Pixi) and ^oxt is an externally 
applied flux. This constraint is used to eliminate 4>2 from 
Eq. jSH). 

In obtaining the Hamiltonian, we assume that the 
qubit does not significantly perturb the resonator such 
that the mode decomposition for ^{x,t) found in the pre- 
vious section is a good approximation even in the pres- 
ence of the qubit. This approximation is accurate for 
small qubit capacitances such that the capacitivc terms 
in Eq. (|3.ip do not induce large frequency shifts of the 
resonator, and if the inductance of the center line of 
the resonator of length w = X2 ~ xi where the qubit 
is connected is smaller than the total inductance of the 
qubit [i.e. L°{xi)w/Yl,k^3:k < 1 with Lj^k = ^o/Ejk 
the Josephson inductance of junction k) such that most 
of the current is flowing through the resonator. Both of 
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these assumptions can safely be satisfied in practice with 
smah junctions. We note that while these constraints 
are useful in deriving the system Hamiltonian, the main 
results can hold even if they are not strictly respected. 

It is useful to introduce the sum and difference fluxes 
<i)± = {[</'3 + ^(2:2)] ± + ip{xi)]}/2, where (j>i + ip{xi) 
and (/)3 + '>p{x2) represent the flux on the island of the 
qubit separated from the resonator by junction 1 and 
3, respectively. The charges conjugate to these fluxes 
are q± = 9£/9(/)±. Using these conjugate variables, the 
Hamiltonian is easily obtained in the usual way [ll| . Af- 
ter transformation under the unitary T+T_, with 



T.=n 



cxp 



(3.3) 



where 5^ = m„(x2) ± u„(xi) and using phase variables 
= 4>tJ./'^0i the Hamiltonian reads 



9: 



2Cn 2C+ d 

Ej [2 cos iy9_|_ cos iy9_ + a C0s((y5ext + <^ + 2(y5_ ) 



In this expression, 1^ is a quantum flux bias given by 

ip^'^tpnS~/(po, (3.5) 

n 

where in general \ip\ <^ 1, and the resonator mode- 
dependent capacitances are 



1 



2Cr + 2C,2iS-y 



Cn CI 

1 _ a(fi-i + 2c^) + QiC,-2(^-)^ 



(3.6) 



with Cl = 2[Cr{Cji + 2Cj2) + {KifCjiC.j2]. The Joseph- 
son potential energy takes the usual form for a 3-junction 
flux qubit jT^ . The usefulness of the unitary transfor- 
mation is to change the phase bias from vacuum fluc- 
tuations in the resonator field to a flux bias <^ directly 
on the qubit which is simply adding to the external flux 

Dcflning the qubit capacitances = 'Ylm ^^'^ 
panding the term proportional to aEj to first order in 
If, the resulting Hamiltonian becomes 



H = Hr- 



H, 



qb 



(3.7) 



where = 'Ylm ^^n^AiO-n is the resonator Hamiltonian 

Hqb 



2C- 2C+ (3.8) 

Ej [2 cos ifj^ cos if^+a cos(iy9ext + 2iy9_ )] 



the standard fiux qubit Hamiltonian and 



Hq = -Y.^S-q^t 



aEjipsin{ipe^t + 2ip^) 



(3.9) 
(3.10) 



describe charge and fiux coupling of the qubit to mode n 
of the resonator, respectively. 

Projecting on the eigenstates {|fc)} of frequencies {ilk} 
of the qubit Hamiltonian Hqb, the fiux coupling Hamil- 
tonian can be expressed as 



n k,l 



where 



%0 „ a£^jA</?„(fc| sin((/?e:,t 4- 2(^_)|/). 



(3.11) 



(3.12) 



Here we have used Eq. 



and defined A(^„ = 



5^ \Jh/2C.rL0n/ '■Po- These matrix elements are easily eval- 
uated after diagonalizing Hqb numerically to find the ex- 
act qubit eigenstates. At the fiux sweet-spot, $cxt — 
(3.4) *I'o/2, only off-diagonal couphng g'^^ between states k 
and / of different parity remain [3]. 

Using the expression of Eq. (|2.9p for the above 
selection rule reduces the charge coupling Hq to 



^4 = E E H\nm{i\ - \i){k\) (at _ a„: 

n k,l>k 

for states |Z) of different parity and where 



(3.13) 



yq,n y 2 



(3.14) 



is a real quantity and maximal between states (fc, I) 
(1, 2). Comparing gq to g^p we get 



q¥ 



{2a + l)Ej {k\ Sin(^e:rt + 2ip_)\l) 



huJr, 



\q-/2e\l) 



(3.15) 



Since in practice Ej ^ fkUn for fiux qubit, the charge 
matrix elements are at best a fraction of unity in the 

vicinity of flux degeneracy point, we find 



k,l / kA 

9a 19, 



> 



10 — 10 . Unsurprisingly, charge coupling is negligible. 



IV. JAYNES-CUMMINGS HAMILTONIAN 



In the rotating-wave approximation (valid when 
<C {uJn,i^min{k,i)})j the fuU Hamiltonian takes the 
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Figure 4: (Color online) Cut along the ip+ = axis of the 
double well potential of the flux qubit with the first (full blue) , 
second (dashed red) and third (dotted green) eigenstates for 
qubit parameters Eji = 259 GHz, a = 0.8 and Ej/Ec ~ 35 
with $cxt = "l>o/2. 



Jaynes-Cummings form |13| 



■E 

n.k.l 



where 5^' = ^r^^ - g-.^^ 



^kl 



%n(|fc)ai4+h-C.) 

(4.1) 



In addition to the strong coupling and the low charge 
noise, an important advantage of studying the Jaynes- 
Cummings physics in this system is the very large an- 
harmonicity of the flux qubit compared to the trans- 
mon 0]. This is illustrated in Fig. 2] which shows the 
first few eigenenergies of the flux qubit Hamiltonian Hg. 
Moreover, with flux qubits it is also possible to take ad- 
vantage of the fact that the first two eigenstates can be 
localized in the wells of the potential, while the higher 
eigenstate is delocalized. These three states can be used 
as a A system Q , opening possibilities for many quantum 
optics phenomena in cavity QED with superconducting 
circuits. 

It is also worth pointing out that multiple qubits can 
be coupled to the same resonator. With the qubits fabri- 
cated in close proximity to a node of the eigenfunction to 
which they are (most strongly) coupled, the Hamiltonian 
of the system simply reduces to 



H = Hr 



N 



(4.2) 



where the coupling Hamiltonians are calculated by pro- 
jecting the operators onto each qubit subspace. Two- 
qubit gates can be generated in this system in the same 
way as with charge-qubit based circuit QED 14 1 . 



V. COMPARISON TO GEOMETRIC 
COUPLING AND NUMERICAL RESULTS 



As shown in the last section, qubit-resonator flux cou- 
pling is provided by the vacuum fluctuations of the res- 
onator field (p threading the qubit loop. For a qubit shar- 
ing a length w with the resonator, this flux can be ex- 
pressed as 



ip = 



E 



wLo{xi)i{xi)/(po 



„(.Ti +w) - u„{xi)] {al + an)/(po 

(5.1) 



where Lo(a:i) and/(a;i) = —dxtpix) / Lq(x)\x-^ are, respec- 
tively, the inductance per unit length and the current in 
the resonator at the location xi of the qubit. The slope 
of the flux field dxipix) at the location of the qubit is 
given in Table [T] for different resonator geometries and 
materials. As can be seen there, coupling is enhanced by 
locally increasing the resonator inductance. This is to be 
compared to the case where the qubit loop is mutually 
coupled to the resonator only by mutual inductance 
In this situation, one finds ip = MI{xi) / ip^ with 



M ' 



2tt 



In 



L 



(5.2) 



as given from the Neumann formulae for a rectangular 
loop of length w and width L separated by a distance d 
from the resonator center conductor, here approximated 
by a infinite cylindrical wire. Geometric coupling will 
win over direct coupling only if the ratio L/d is made 
such that 



— > e fo/^T 
d ^ 



1. 



(5.3) 



Since L^{xi) can reach several units of /io/27r in the con- 
striction due to the contribution of the kinetic induc- 
tance, geometric coupling can only win by either increas- 
ing the qubit area (w x L) or by reducing the distance 
between the qubit loop and the resonator central line (d) . 
Large loops will make the qubit more susceptible to sur- 
rounding flux noise, while placing the qubit very close to 
the resonator can be challenging in addition to increasing 
the capacitive coupling to the resonator. In contrast, di- 
rect coupling leads to a coupling strength that scales with 
length w of the shared part between qubit and resonator 
rather than with the area. Strong coupling can therefore 
be reached without large sensitivity to flux noise. 

Table H] shows the coupling g^^^ of directly connected 
qubits for various transmission-line configurations. To 
obtain these results, we have taken parameters close to 
those of Ref. [H with Eji = 259 GHz, a = 0.8 and 
Ej/Ec = 35. These parameters where used in Fig. 2] 
The qubit 's loop width was taken to he w = 5 fim. 
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1000 3000 5000 7000 9000 

Ej (GHz) 

Figure 5: Flux bias Atp\ and resultant qubit-resonator cou- 
pling energy g^^-^ as calculated from Eq. (|3.12p induced by a 
Josephson junction placed in the center of a homogenous alu- 
minum transmission-line resonator (see Table |lj as a function 
of the Josephson energy Ej. The inset shows that the qubit- 
resonator coupling energy can reach several tens of percent of 
the resonator frequency. Qubit parameters are given in the 
text. 



For these realistic values, this system easily reaches the 
strong coupling regime. Moreover, in all cases shown 
here, the coupling strength is substantially larger than 
the estimates for geometric coupling d, Q . 



VI. ULTRA-STRONG COUPLING REGIME 

Using a Josephson junction to locally change the in- 
ductance can result in much stronger coupling. As il- 
lustrated in Fig. [5l for a relatively large Josephson ener- 
gies Ej > 1000 GHz, the phase bias seen by the qubit 
is so large that the coupling energy can easily reach 
g^^^/2Ti ^ 1000 MHz and beyond, corresponding to sev- 
eral tens of percent of the resonator frequency. This cou- 
pling can be increased further by lowering the Josephson 
energy of the inserted junction, as long as the correspond- 
ing Josephson inductance is small compare to that of the 
qubit. 

In this ultra-strong coupling regime 0, [i3j the RWA, 
used in going from Eq. (|3.1ip to Eq. (|4.ip . breaks down 
and the full Hamiltonian must be considered. This cir- 
cuit then becomes a solid-state implementation of the 
E ® (3 Jahn- Teller model . While coupling of the 
artificial atom to the electric field of the resonator docs 
not lead to super-radiant phase transition, magnetic cou- 
pling, which dominates here, does [H, [l^. In this situ- 
ation, the ground state of the combined qubit-resonator 
system can be an entangled state corresponding to the os- 
cillator being displaced by a qubit-state dependent quan- 
tity [21], |2^. In this ground state, we can expect a finite 
photon population, something which could be measured 
using a second resonator in a number splitting experi- 
ment ii,!!!. 



In addition to the breakdown of the RWA, higher- 
order terms in the expansion of the Josephson energy in 
Eq. (|3.4p have to be taken into account as the coupling 
strength increases. Second-order corrections lead to an 
additional flux coupling Hamiltonian H ~ of the form 



(6.1) 



where 



kl 



-EjiAipn)^{k\ cos{^ext + 2^-)|0 (6.2) 



This second-order correction leads to ac-Stark shifts and, 
more interestingly, can be used to generate squeezing of 
the microwave field inside the resonator. Tuning the 
circuit parameters can lead to detectable effects with 
C|!„/27r > 1 MHz. We note that this also leads to 
resonator mode-mode coupling. In practice however, 
the frequency separation between these modes is large 
enough that this can be neglected. 



VII. CONCLUSION 

We have obtained the Hamiltonian of a superconduct- 
ing flux-qubit directly coupled to the center conductor 
of a coplanar transmission-line resonator. By using a 
constriction in the center line of the resonator, the cou- 
pling strength between the qubit and the resonator can 
be significantly increased. This is due to the increase 
in the geometric and kinetic inductance of the line and 
the resulting large phase bias seen by the qubit. There 
arc two main advantages of this approach compared to 
coupling based on the mutual inductance between the 
qubit and the resonator: the coupling is much stronger 
in magnitude and this is possible without working with 
large qubit loops which would increase sensitivity to flux 
noise. Together with the insensitivity of flux qubits to 
charge noise, its large anharmonicity and its A config- 
uration, this approach leads to the possibility of study- 
ing numerous quantum optics effects with superconduct- 
ing circuits. Finally, by replacing the constriction with 
a Josephson junction of large Josephson inductance, we 
have shown that the coupling can be as large as several 
tens of percent of the resonator frequency. In this sit- 
uation, the breakdown of the RWA should lead to an 
entangled qubit-resonator ground state. 
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a) 



S(x) 









b) 



c) 




the central line \2l 



Figure 6: Geometry of the inhomogeneous transmission line 
used in numerical simulations. Central-line width S{x) and 
ground- plane spacing W{x) are shown in a) for the total 
length of the line 21 — h mm, and in b) around the con- 
striction indicated by the dashed square in a). The thickness 
of the central-line is shown in c). In this example, the outer 
dimensions of the homogeneous part of the line are Sout — 5 
^tm, WoMt = 2.5 pim and tout = 200 nm while the inner di- 
mensions inside the constriction at a: = with 5in = 50 nm. 
Win = 4.95 ^ra and tin = 50 nm. 



Appendix A: DESIGN AND CHARACTERISTICS 
OF INHOMOGENEOUS TRANSMISSION LINES 



Appropriate modeling of the inhomogeneous 
transmission-line electrical characteristics is needed 
to compute eigenmodes, frequencies and ultimately 
the coupling between the qubit and the resonator. In 
this section, we give details on the geometry of the 
inhomogeneous resonators that were used in Table [H 

The capacitance per unit length C*' {x) , the inductance 
per unit length L^{x) and the impedance Z°{x) of the 
coplanar transmission-line resonator depend on the ratio 
between the width of the center electrode S and the dis- 
tance between the two ground planes S + 2W, W being 
the distance between the ground plane and the edge of 



goo 



-'^kin 



(Al) 



CO 



where is the dielectric constant of t he subt rate, fco = 
S/ [S + 2W) is the aspect ratio, fcg = ^^1 — and K{x) 
is the complete elliptic integral of the first kind. By de- 
creasing the aspect ratio fcg along the line, the induc- 
tance and impedance of the line are locally increased 
while the capacitance is decreased. For a superconduct- 
ing resonators, the kinetic inductance can be expressed 

as m 



^kin =MoAl(T) 



C 



4:ADK{ko) 
0.4 



1.7 



sinh(i/2AL(r)) 



^[iB/Ar-l][l-iB/Dy 



(A2) 



where Xl{T) is the London penetration depth of the su- 
perconductor at temperature T, t is the thickness, and 



A: 




B 



W^, D = 



2t 



(A3) 



C. 



As for the geometrical inductance, a decrease of the as- 
pect ratio will increase the kinetic inductance of the line 
but the effect is rather marginal unless the dimensions of 
the cross section of the central electrode become of the 
order of A^. 

The geometry of the inhomogeneous transmission-line 
used in Table U is depicted in Fig. We consider 

a regular, initially homogeneous transmission line res- 
onator made of cither aluminum (Al = 16 nm) or nio- 
bium (Al = 39 nm) with a total length of 2£ = 5 mm. 
The central electrode is t = 200 nm thick and S = 5 ^m 
wide. The ground planes are W = 2.5 /im away from the 
edge of the central-line. For simplicity, the distance be- 
tween the ground planes {S + 2W) is not modified. The 
width and thickness of the central electrode are reduced 
at the center of the resonator to create a constriction. 
The dimensions are reduced gradually over a length d 
from the initial outer values S'out , Wont and tout down to 
the minimal inner values 5*111, and tin at a; = 0. The 
electrode width S{x), thickness t{x) and ground-planes 
spacing W{x) are continuous smooth functions for ease 
of computation. 

Table U summarizes the numerical results obtained for 
inhomogeneous resonators made of Al and Nb and char- 
acterized by width, thickness and ground-plane spac- 
ing as illustrated in Fig. [Sj For aluminum, as the 
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Figure 7: (Color online) First mode (normalized to 1) for the 
niobium transmission line (dashed black) and for a homoge- 
neous aluminum transmission-line resonator intersected by a 
(pointlike) Josephson junction of plasma frequency lo^jl-K — 
40 GHz with Ej = 6000 GHz (full blue). The transmission- 
line characteristics are detailed in Table U] The inset shows an 
enlarged view of the mode in the vicinity of the junction. The 
abrupt variation in the mode at the position of the junction 
enables a much greater flux bias on the qubit. 



central line cross section dimensions are reduced from 
5 /im X 200 nm down to 50 x 50 nm^, the total in- 
ductance per unit length can be increased by a factor 
of 4. As it shown in Fig. [21 the slope of the flux field 
|92;V'i(a;)| = \dxUi (x) \ \J h/2CrOJi inside the constriction 
is also increased by a factor of 4. On the other hand, be- 
cause of the larger London penetration depth, the kinetic 
inductance of niobium resonators can be very important. 
In this case, the local inductance can be increased by a 
factor of '-^ 8 for a cross-section of 50 x 50 nm^, leading 
to a fivefold increase in the slope of the fiux field. 



description of a resonator with an integrated Josephson 
junction can be found elsewhere [loj and we recall only 
the main results (see Fig. [7|). 

For a transmission-line resonator interrupted by a 
Josephson junction of linear Josephson inductance Lj = 
(Pq/Ej and capacitance Cj, the cigcnmodcs can be de- 
scribed by spatially oscillating functions Um{x) given by 



cos[fc,n(a; + i)] X < Xj 

Bm COs[k,n{x — i)] X > Xj 



(Bl) 



If the junction is placed at the center a; = of the res- 
onator, Brn = — 1 and the wave- vectors km are solutions 
of the transcendental eigenvalue equation 



^(1-^1 cot[fc,„^] 



(B2) 



where = fcm/ v L^C^ are the resonance frequencies of 
the circuit and ujp = 1/^/LjCj is the plasma frequency 
of the junction. The eigenmodes of the circuit are found 
to obey a generalized orthogonality equation 



dxCoUjnix)ujn' (x) + Cj5^5^^^, = Ct,5„ 



(B3) 



where iS" — Wm(0"'") — Um(0~) is the dimensionless mode 
gap at the junction and = 2£C° -f Cj is the total 
capacitance of the circuit. This orthogonality equation 
is used to fix the normalization A^- 

It follows that in the linear approximation of the 
Josephson inductance, the resonator can be describe by 
sum of harmonic oscillators of frequency uom 



{alam + 1/2) 



(B4) 



Appendix B: TRANSMISSION LINE 
INTERSECTED BY A JOSEPHSON JUNCTION 



Inserting a Josephson junction in the center conduc- 
tor of the resonator at the location of the qubit can lead 
to significantly stronger coupling. We note that having 
a four, rather than three, junctions is natural for a flux 
qubit [23| . If the Josephson inductance of the resonator 
junction is much smaller than the total inductance of the 
qubit loop, the qubit acts once again as a simple per- 
turbation on the resonator cigcnmodcs. The theoretical 



where the ladder operators obey the commutation rela- 



tion [a,i 



and define fiux ijjm and charge q„ 



operators given by Eq. (|2.9p with the appropriate capaci- 
tance and frequency definitions. As it is shown in Fig. [3 
the presence of the junction creates a very abrupt dis- 
continuity in the modes. This leads to a large fiux field 
slope dxtpi{x), which in turns leads to very strong qubit- 
resonator coupling. Wc note that very strong coupling 
can be attained even for negligible non-linearity of the 
resonator mode (l0| . 
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